arXiv: 1504.04269vl [quant-ph] 16 Apr 2015 


Fate of Accidental Symmetries of the Relativistic 
Hydrogen Atom in a Spherical Cavity 


M. H. Al-Hashimi“’^, A. M. Shalaby“’^, and U.-J. Wiese^*^ 

“ Department of Mathematics, Statistics, and Physics 
Qatar University, A1 Tarfa, Doha 2713, Qatar 
^ Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics 
Bern University, Sidlerstrasse 5, CH-3012 Bern, Switzerland 
Physics Department, Faculty of Science, Mansoura University, Egypt 

INT, Washington University, Seattle, U.S.A., Preprint: INT-PUB-15-009 

PACS numbers: 


April 17, 2015 

Abstract 

The non-relativistic hydrogen atom enjoys an accidental S'0(4) sym¬ 
metry, that enlarges the rotational S'0(3) symmetry, by extending the 
angular momentum algebra with the Runge-Lenz vector. In the relativis¬ 
tic hydrogen atom the accidental symmetry is partially lifted. Due to 
the Johnson-Lippmann operator, which commutes with the Dirac Hamil¬ 
tonian, some degeneracy remains. When the non-relativistic hydrogen 
atom is put in a spherical cavity of radius R with perfectly reflecting 
Robin boundary conditions, characterized by a self-adjoint extension pa¬ 
rameter 7 , in general the accidental S'0(4) symmetry is lifted. However, 
for R = (/ -|-1)(/ -|- 2)a (where a is the Bohr radius and I is the orbital an¬ 
gular momentum) some degeneracy remains when 7 = cx) or 7 = ^. In 
the relativistic case, we consider the most general spherically and parity 
invariant boundary condition, which is characterized by a self-adjoint 
extension parameter. In this case, the remnant accidental symmetry 
is always lifted in a hnite volume. We also investigate the accidental 
symmetry in the context of the Pauli equation, which sheds light on 
the proper non-relativistic treatment including spin. In that case, again 
some degeneracy remains for specihc values of R and 7 . 
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1 Introduction 


Already in 1873, Bertrand has shown that the 1/r and potentials are the only 
spherically symmetric potentials for which all bound classical orbits are closed pQ. 
Accidental dynamical symmetries arise in many quantum mechanics problems [2], 
ranging from the harmonic oscillator to the hydrogen atom and to a charged particle 
in a constant magnetic field. In 1935 Fock noted that the hydrogen atom possesses 
a “hyper-spherical” S'0(4) symmetry [3], and Bargmann |1] showed that the gener¬ 
ators of the accidental symmetry are the components of the Runge-Lenz vector [5]. 
For an isotropic harmonic oscillator in d dimensions, the spatial rotational symmetry 
SO{d) is dynamically enhanced to SU{d) and for the non-relativistic hydrogen atom 
it is enhanced to SO{d -|- 1). For a charged particle in a constant magnetic held, 
the center of the circular cyclotron motion plays the role of the Runge-Lenz vector 
In this case, translation invariance (up to gauge transformations) disguises 
itself as an accidental symmetry. In all these cases, at the classical level the acciden¬ 
tal symmetry implies that all bound classical orbits are closed curves, while at the 
quantum level it leads to enlarged degeneracies of the energy levels of bound states. 
For a charged particle in a constant magnetic held, the corresponding Landau levels 
are even inhnitely degenerate. 

In the past, we have investigated how accidental symmetries are ahected when 
the geometry of the problem is modihed. For example, in [9j we have studied the 
harmonic oscillator and the hydrogen atom on a 2-d cone with dehcit angle 6. For 
general values of S, the accidental SU{2) = SO{3) symmetry is then completely 
lifted. However, if ^ is a rational number, some accidental degeneracies remain. Re¬ 
markably, the degeneracies correspond to fractional (neither integer nor half-integer) 
“spin” of the accidental SU (2) symmetries. This subtle ehect arises because on the 
cone the Runge-Lenz vector is no longer self-adjoint in the domain of the Hamil¬ 
tonian. Hence, accidental symmetry operations induced by the Runge-Lenz vector 
in general lead out of the domain of the Hamiltonian, and thus no longer repre¬ 
sent physical symmetries. For rational deficit angles, on the other hand, repeated 
applications of the Runge-Lenz vector eventually lead back into the domain of the 
Hamiltonian and thus give rise to remnant degeneracies. In [H] we have changed the 
geometry of the Landau level problem by confining the charged particle to a 2-d spa¬ 
tial torus, threaded by n$ units of quantized magnetic flux. The spatial boundary 
conditions are then characterized by two self-adjoint extension parameters (related 
to the two holonomies of the torus) which explicitly break translation invariance to 
a discrete Z(n#) symmetry. The continuous accidental symmetry of the problem in 
the inhnite volume is then reduced to a hnite magnetic translation group and the 
inhnite degeneracy of the corresponding Landau levels is reduced to n$. 

Spatial boundary conditions are naturally characterized by families of self-adjoint 
extension parameters mm For example, the general Robin boundary conditions 
of a non-relativistic quantum mechanical particle in a box with perfectly reflecting 
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walls are characterized by a real-valued self-adjoint extension parameter 7 , which 
interpolates continuously between Dirichlet and Neumann boundary conditions [Hh 
[13]. In [13] we have used the theory of self-adjoint extensions to derive a generalized 
Heisenberg uncertainty relation for hnite-volume quantum dots with general bound¬ 
ary conditions. Perfectly reflecting boundary conditions for relativistic fermions 
described by the Dirac equation have also been investigated in na. They are char¬ 
acterized by a 4-parameter family of self-adjoint extensions. The same is true for the 
Pauli equation, that arises from the Dirac equation in the non-relativistic limit. In 
[13] we have investigated the supersymmetric descendants of self-adjointly extended 
quantum mechanical Hamiltonians. In particular, we found that the superpartners 
of Hamiltonians obeying general Robin boundary conditions characterized by the 
parameter 7 always obey simple Dirichlet boundary conditions, while 7 determines 
the value of the superpotential at the boundary. 

Hydrogen atoms conhned to a hnite volume have been investigated to mimic the 
effects of high pressure, in order to better understand, for example, white dwarf 
stars. In order to model hydrogen at high pressure, a hydrogen atom at the center 
of a spherical cavity has been studied by Michels, de Boer, and Bijl as early as 1937 
ini- This work was extended by Sommerfeld and Welker in 1938 [18] as well as 
in [Tn]j22]- The effect of general Robin boundary conditions (characterized by the 
real-valued self-adjoint extension parameter 7 ) on the accidental symmetry of the 
non-relativistic hydrogen atom and the harmonic oscillator conhned to a spherical 
cavity have been investigated in [23l [23]. In both cases, in general the accidental 
symmetry is lifted. Again, this is because an application of the Runge-Lenz vector 
leads out of the domain of the Hamiltonian. However, as for the corresponding 
systems on the cone, repeated application of the Runge-Lenz vector leads back into 
the domain of the Hamiltonian, which gives rise to a hnite-volume remnant of the 
accidental symmetry, at least for particular radii of the spherical cavity and for 
specihc values of the self-adjoint extension parameter [231I25H27]. 

In this paper, we extend this analysis to the relativistic hydrogen atom de¬ 
scribed by the Dirac equation [2HH31]- The accidental symmetry of the non- 

relativistic problem is then substantially reduced, but not lifted completely. This 
is due to the Johnson-Lippman operator [3ll [35l [36] , the relativistic counter part 
of the Runge-Lenz vector, which still leads to enhanced accidental degeneracies. 
When the relativistic hydrogen atom is placed inside a conhning spherical cavity, 
the accidental symmetry is completely lifted, because an application of the Johnson- 
Lippman operator again leads outside the domain of the Hamiltonian. In contrast to 
the non-relativistic case, repeated applications of the Johnson-Lippman operator do 
not lead back into the domain of the Hamiltonian, and thus no remnant accidental 
symmetry persists in a hnite volume. The non-relativistic limit of the Dirac hydro¬ 
gen atom leads to the Pauli equation (and not to the Schrodinger equation without 
spin). Interestingly, the conserved current induced by the Dirac equation contains a 
spin contribution. This again gives rise to self-adjoint extensions, which in general 
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induce spin-orbit couplings via the spatial boundary conditions. Remarkably, in this 
case again a hnite-volume remnant of the accidental symmetry persists for particular 
radii of the conhning cavity and for particular values of the self-adjoint extension 
parameter. Our study further illuminates the subtle effects of self-adjoint extension 
parameters on accidental symmetries, in particular, for relativistic quantum systems 
and their non-relativistic counter parts. 

The rest of the paper is organized as follows. Section 2 reviews the accidental 
symmetries of the hydrogen atom, both in the non-relativistic Schrodinger and in 
the relativistic Dirac treatment. In Section 3 we discuss the self-adjoint extension 
parameters characterizing perfectly reflecting cavity walls both in the relativistic 
case and in the non-relativistic case with and without spin. In Section 4 we place 
the hydrogen atom in a spherical cavity. After reviewing the results of the non- 
relativistic Schrodinger treatment, we proceed to the Dirac equation and its non- 
relativistic Pauli equation limit. Finally, Section 5 contains our conclusions. 


2 Accidental Symmetries of the Hydrogen Atom 


In order to make the paper self-contained, in this section we review the accidental 
symmetries of both the relativistic and the non-relativistic hydrogen atom. In the 
non-relativistic Schrodinger problem the accidental symmetry is generated by the 
Runge-Lenz vector, while in the relativistic Dirac problem it is generated by the 
Johnson-Lippmann operator. 


2.1 Accidental Symmetry of the Non-relativistic Schrodinger 
Atom 


Let us consider the Schrodinger Hamiltonian for the hydrogen atom 


H 



( 2 . 1 ) 


Here M is the electron mass and e is the fundamental electric charge unit. Through¬ 
out this paper we put h = 1 but leave the velocity of light c explicit in order to 
facilitate considerations of the non-relativistic limit. Thanks to the isotropy of the 
Coulomb potential, the Hamiltonian obviously commutes with the angular momen¬ 
tum L = r X p, i.e. [H, L] = 0. It is well known, but not so obvious, that the 
Hamiltonian also commutes with the Runge-Lenz vector 


R = 


2M 


p X L — L X p) — e^Cr, 


Cr — “ • 


( 2 . 2 ) 
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The angular momentum L and the Runge-Lenz vector R generate an S'0(4) exten¬ 
sion of the SO{3) rotational symmetry, with the following commutation relations 

‘2H 

[Tj, Tj] [Z/j, Rj] iEijkRfc, i ^ Eij^L^. (2.3) 


The two Casimir operators of the S'0(4) algebra are 


Cl 



Me^ 


C2 


L-R. 


(2.4) 


Using the explicit form of R one can show that L-R = 0, which implies that 
only particular representations of S'0(4) can be realized in the hydrogen atom. 
These representations are n^-fold degenerate, which indeed represents the correct 
accidental degeneracy of the hydrogen atom bound states. For example, the 2s state 
is degenerate with the three 2p states, and the 3s state is degenerate with the three 
3p and the hve 3d states. The value of the Casimir operator in the corresponding 
representations is Ci = n^, such that 


H 


Me^ ^ Me^ 

- ^ = - 

2Ci 2n2 ’ 


(2.5) 


which indeed is the familiar spectrum of hydrogen atom bound states. 


2.2 Accidental Symmetry of the Relativistic Dirac Atom 


Let us now proceed to the Dirac Hamiltonian 

e2 

H = a ■ pc-\- (3Mc^ -, 

r 

using the following conventions for the Dirac matrices 




( 2 . 6 ) 


(2.7) 


Here 1 and 0 represent 2x2 unit- and zero-matrices, and a are the Pauli matrices. 
The Hamiltonian commutes with the total angular momentum 

j = L + S = L + ^, [H,J]=0. (2.8) 

Furthermore, it is parity invariant and thus it commutes with the operator 

P = /3J, [H,P] = 0, [J,P]=0, (2.9) 


where / performs a spatial inversion. Dirac has also discovered another symmetry 
of the system generated by the operator 

K = /3{L-L + 1), [H,K]=0, [J,K]=0, [P,K] = 0. (2.10) 
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The eigenvalues k of K follow from 


a-L + l = {L + Sf-L^-S^ + l = P-L^ + ^ ^ 
k = j{j + 1) - /(/ + 1) + ^ = j{j + 1) - (^j T 0 (•^' ^ ^ + ^) + i 

= ±(j + ^). (2.11) 

The sign of the quantum number k indicates whether L and S couple to j = / + ^ 
{k = i + \) OT to j = I — \ {k = —(j + ^)). Energy eigenstates are characterized by 
the principal quantum number n as well as by j, ja, and fc, such that 

H\njhk) = Enk\niji.k), J = j(j + l)|njj3/c), .h\nji^k) = 

KWii^k) = k\nii:ik). (2.12) 


Note that in some other works the eigenvalue of K is denoted as —k. The energy 
eigenvalues (of positive energy states) are given by [2H1 [29] 


Enk — Me ( 1 + 


a 


(n — \k\ + 


- 1/2 


(2.13) 


Here a = ^ jc is the hne-structure constant (in units were h = 1). In the non- 
relativistic limit one recovers the Schrodinger result, corrected by spin-orbit hne- 
structure effects 


Enk = Mc^ 


2n^c^ 


e® / 1 

2n®c^ \j + I 



(2.14) 


The hne-structure leads to a partial lifting of the large n^-fold degeneracy related 
to the accidental S'(9(4) symmetry of the non-relativistic problem (without spin). 
In particular, the energy no longer just depends on n but also on j and, for a given 
n, is slightly larger for the states with larger j. Still, a remnant of the accidental 
symmetry persists even for the Dirac equation, because the energy depends only on 
\k\. As a result, the states with j = I + ^ {k = j + ^), which have orbital angular 
momentum I, and the states with j = I + 1 — ^ [k = —(j -|- 1)), which have orbital 
angular momentum Z -|- 1, are degenerate (as long as they have the same principal 
quantum number n). This gives rise to an accidental 2(2j -|- l)-fold degeneracy. The 
states with the maximal orbital angular momentum I = n — 1 and the maximum total 
angular momentum j = l + \= n — only have the usual (2j -|- l)-fold degeneracy 
that results from the S'0(3) rotational symmetry. For example, the 4 states 2Si/2 
and 2 Pi /2 with total angular momentum j = \ (in the standard spectroscopic nlj 
notation) are degenerate while the 4 states 2 P 3/2 with j = \ have a slightly higher 
energy. Similarly, the 4 states 3Si/2 and 3Pi/2 are accidentally degenerate, the 8 
states 3 P 3/2 and 3 D 3/2 are again accidentally degenerate, but the 6 states 3 D 5/2 with 
the maximal J = 3 ~ | § are not. 
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The relativistic counterpart of the Runge-Lenz vector, which again commutes 
with the Hamiltonian, is the pseudo-scalar Johnson-Lippmann operator [35l [36] 


A = iK'j^ 


H 


{^Mc^ 


13 ] — aS • Cr, [-H, A\ = 0, [J, A\ = 0. 


(2.16) 


The operator A anti-commutes with K, i.e. {K,A\ = 0, while A^, which plays the 
role of a supersymmetric “Hamiltonian” PTIHO] is directly related to H and K via 







“t“ (y. 


(2.16) 


The eigenvalues of A^ are hence given by 


a" = 






I 2 2 

(y — (y — 


\Mc^ J 

The operator A acts on energy eigenstates as 

Alnjjik) = a|njj3 - k), 


{a — | 1 :| + y/k^ — + a'^ 


(2.17) 


(2.18) 


i.e. it relates two accidentally degenerate states with quantum numbers 3ik. If a 
state has maximal j = n — it is annihilated by the operator A because then 
n = j + ^ = k=\k\ such that 

- —---- = 0 ^ a = 0. (2.19) 

k^ — 

As was pointed out in |39|, the Johnson-Lippman operator generates a supersymme¬ 
try that relates the accidentally degenerate states. The unpaired state with maximal 
j = n — I is the lowest state underneath a tower of paired states with the same value 
of j but higher values of n. The fact that the lowest state is not paired indicates 
that supersymmetry is not spontaneously broken. 


3 Perfectly Reflecting Cavity Boundary Condi¬ 
tions 


Following na, in this section we discuss perfectly reflecting boundary conditions, 
first in the context of the non-relativistic Schrodinger equation and then for the 
relativistic Dirac equation. Next we take the non-relativistic limit of the Dirac 
equation and arrive at the Pauli equation. The spin then enters the conserved 
current, with important implications for the self-adjoint extension parameters that 
characterize the spatial boundary conditions. 
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3.1 Non-relativistic Robin Boundary Conditions 


Let us consider the general Schrodinger Hamiltonian 


H 


2M 


+ V{x) 


1 


A + V{x). 


(3.1) 


The continuity equation that guarantees probability conservation then takes the 
form 

dtp{x,t)+ V ■ j{x,t) = 0, (3.2) 

where the probability density and the corresponding current density are given by 


p{x,t) = \^{x,t)\^, j{x,t) 


1 

2Mi 


T(x, t)*V\h(a;, t) — VT(T, t)*T(a;, t) 


(3.3) 


We now consider an arbitrarily shaped spatial region H and we demand that no 
probability leaks outside this region. This is ensured when the component of the 
probability current normal to the surface vanishes 

n{x) ■ j{x) =0, xG 911. (3.4) 

Here n{x) is the unit-vector normal to the surface dQ at the point x. The most 
general so-called Robin boundary condition that ensures probability conservation is 
given by 

7 (x)T(x)-I- fi(x) ■ VT(x) = 0, X e dfl, (3.5) 

which indeed implies 

n{x) ■ j{x) = 4 '(x)*7(x)T(x, t) 4- 7 (x)*T(x, f)*\k(x, t)] =0, x G 911, (3.6) 

such that 7 (x) G M. 


It is important to show that the Hamiltonian endowed with the boundary con- 

se: 

{X\H\'^) = Id^xxixy 


dition eq.(3.5) is indeed self-adjoint. First, we investigate 

1 


in 


2M 


A + V{x) 


T(x) 


(TX 


in 


1 


2M 


Vx(x)* ■ VT(x) -1- x(x)*H(x)\['(x) 


2M „ 

f d^x 

Jn 


an 


dn ■ x(x)*VT(x) 


2M 


A + V{x) 


X{xy \ T(x) 


2M 


dn ■ 


'an 


= mH\xr + 


Vx(x)*\k(x) — x(^)*V\k(x) 
1 


2M 


dn ■ 


Ian 


Vx(x)*T(x) — x(x)*VT(x) . (3.7) 

























Consequently, the Hamiltonian is Hermitean (or symmetric in mathematical par¬ 
lance) if 


dn ■ 'Vx{x)*'^{x) — x{x)*'V'^{x) =0. 
an L 

Using the boundary condition eq.( |3.5 ), eq.( |3.8 ) simplihes to 

/ d'^x n{x) ■'Vxix)* +'yi^)xid^)* '^{x)=0. 
Jan 


(3.8) 


(3.9) 


\h(x) itself is not restricted at the boundary, and hence the Hermiticity of H requires 


n{x) ■ Vx{x) + 7 (x)*x(x) = 0. (3.10) 

For 7 (x) G M, this is again the boundary condition of eq.( |3.5[) . Since x(F) must 
obey the same boundary condition as th(a;), the domain of H\ coincides 

with the domain of if, D{H). Since D{W) = D{H), the Hamiltonian is not only 
Hermitean but, in fact, self-adjoint. 


3.2 Relativistic Cavity Boundary Conditions 

A Dirac particle in a 1-d box has been considered in [H] 02] , and general boundary 
conditions have been investigated in |13]. In [T5| we have investigated general per¬ 
fectly reflecting boundary conditions for Dirac fermions in a 3-d spatial region D, 
thereby generalizing the standard MIT bag boundary conditions |11H16|. 

To keep the discussion as general as possible, we investigate Dirac fermions 
coupled to an external static electromagnetic field, such that the Hamiltonian is 
given by 

H = a ■ (j)c -I- eA(T) j -f- /SMc^ — ed>(T) = —ia ■ Dc + (dMc^ — e^{x). (3.11) 

Here <h(a;) is the scalar and A{x) is the vector potential. The covariant derivative 
then takes the form 

D = W + i-A(x). (3.12) 

c 

The Hamiltonian acts on a 4-component Dirac spinor T(a;, t). In the relativistic 
case, the continuity equation is given by 

dtp{x,t)+ V ■ j{x,t) = 0, (3.13) 

where 

p(T, f) = \k(a;, f)'*'T(a;, f), j{x,t) = c\k(x, f)'*'a\k(a;, f). (3.14) 

Under time-independent gauge transformations, the gauge and fermion fields trans¬ 
form as 

‘^<I)(x) = <F(T), ‘^A{x) = A{x) — Vp{x), ‘^\k(T) = exp (^i-Lp{x)^ 4/(1;). (3.15) 
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Just as in the non-relativistic case, we investigate the Hermiticity of the Hamil¬ 
tonian by considering 

{x\H\'^) = d^x x{x)^ d ■ (—icV + eA{x)) +/3Mc^ — e^{x) \l/(x) 


d^x 


d ■ [—icV + eA{x)) + (3Mc^ — e^{x) x(x) i 'h(x) 


— ic dn ■ x{xYd'^(x) 

Jan 

= {'^\H\xY — ic / (Jn • x(x)'*'aJ/(x). 

JdQ 

This implies the Hermiticity condition 

x(x)^n(x) ■ a\h(x) = 0, x G Shi. 

The corresponding appropriate self-adjoint extension condition is given by 


T3(x) 

T4(x) 


= 


^i(x) 

T2(x) 


A(x) G GL(2,C), xedn, 


(3.16) 

(3.17) 

(3.18) 


which reduces eq.(|3.17) to 

n(x) ■ a 




it 


0 


\l/(x) = 


n[x) ■ a 0 

[(Xi(^)*, X 2 {x)*) n{x) ■ dX{x) + (X 3 (^)*, x4^T) n{x) ■ a] 
Xsix) \ xi(x) 


Ti(x) 

T2(x) 


= —n(x) ■ dXixynix) ■ a { 

XAx) ) \ \ ^ 

Self-adjointness of H requires that D{H) = D{HY, which demands 

A(x) = —n{x) ■ c7A(x)^n(x) ■ a n{x) ■ c7A(x) = — [n(x) ■ aA(x)]^ 


= 0 


(3.19) 


(3.20) 


Consequently, n{x) ■ dX{x) must be anti-Hermitean. This means that, in contrast 
to the non-relativistic Schrodinger problem, relativistic Dirac fermions have a 4- 
parameter family of self-adjoint extensions characterizing a perfectly reflecting wall. 
In the standard MIT bag model [441 - 146] the boundary condition corresponds to 
A(x) = in{x) ■ d. This is a natural choice, because it maintains spatial rotation 
invariance around the normal n(x) on the boundary, but it is not the most general 
possibility. As it should be, the self-adjointness condition eq.(3.18) is gauge covariant 
and it ensures that 


f?(x) ■ j(x) = c\h(x) 


0 n(x) ■ a 


T(x) = 


n{x) ■ a 0 

= c [(\l>i(x)*, \I/ 2 (x)*) n(x) ■ dX{x) + (T 3 (x)*, T 4 (x)*) n{x) ■ a] 


= c (\l/i(x)*, \h 2 (x)*) [n(x) ■ (jA(x) -|- A(x)'*^n(x) ■ a] 


Ti(x) 

T2(x) 


Ti(x) 

T2(x) 


= 0. (3.21) 
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3.3 Non-relativistic Boundary Conditions with Spin 


Following [T5j, we now consider the non-relativistic limit, in which the lower com¬ 
ponents of the Dirac spinor reduce to 


( d'3(f) 

V 



2Mci 


-a ■ D^>{x). 


(3.22) 


Here the 2-component Pauli spinor is given by 


\['(F) = 


V ^2(^) 



(3.23) 


To leading order, the Dirac Hamiltonian then reduces to the Pauli Hamiltonian 

(pc + eA{x)^ 

H = Mc^ + - - — - - e^{x) + pa ■ B{x), (3.24) 

2Mc^ 

with B{x) = V X A{x) being the magnetic held and p = ej^Mc being the Bohr 
magneton, i.e. the magnetic moment of the electron. Here we neglect higher order 
contributions such as spin-orbit couplings and the Darwin term. 


The self-adjoint extension parameters that characterize the most general per¬ 
fectly rehecting boundary condition crucially depend on the form of the conserved 
current. Therefore we now use eq.(3.22) to obtain 


J(f) = c(Ti(f)*,T2(f)*,T3(x)*,T4(f)*) 


0 a 
a 0 


/ ^i(^) \ 

T 2 (f) 

T3(f) 

V ^4(5^) y 


2Mi 

1 

^ L 
1 


T(f)'^DT(x) - (DT(x))'^T(x) 

T(x, t)'*'(T X Z)\I/(x, f) — (Z3\l/(x, t))^ X (jT(x, t) 
^/(x, f)^D\h(x, f) — (Z3\l>(x, t))^\h(x, t) 


2Mi L 

X [T(x,f)^aT(x,t)] . 


(3.25) 


Besides the non-relativistic probability current that is familiar from the Schrodinger 
equation, in the context of the Pauli equation an additional spin contribution arises 
from the Dirac current. Since it is a curl, the spin contribution is automatically 
divergenceless, and the continuity equation is given by 


dtp{x,t) + V ■ j(x,t) = 0, 


(3.26) 
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with the probability density p{x,t) = t). Although the current would 

also be conserved without the spin term, this contribution naturally belongs to the 
current that emerges from the Dirac equation. 

Again following [T^, we now consider the gauge covariant boundary condition 


7 (x)\b(a:)+n(af)- {x) — ia x D'^{x) =0, 'y{x) e GL{2,C), x E dQ, (3.27) 

which implies 


n[x) ■j[x) = 


2Mi 

1 

^ L 
1 

2Mi 


\b(a;)'*'n(x) ■ l3\['(x) — (n(x) ■ l3\I'(x))'^\['(x) 


'^{x,tyn ■ fa X D'^{x,t)] — n ■ f (Dtbfx, t))' x aj 


[—'l/(x)'^7(x)\['(x) + \b(x)''’7(x)'^\['(x)] = 0, 


such that 


7(f)t = ^(f). 


(3.28) 


(3.29) 


As in the fully relativistic Dirac problem, also in the non-relativistic Pauli prob¬ 
lem with spin we again obtain a 4-parameter family of self-adjoint extensions, now 
parametrized by the 2x2 Hermitean matrix 7 (x). 


As in the Schrodinger and Dirac cases, partial integration leads to the Hermiticity 

(3.30) 


condition for the Pauli Hamiltonian of eq.(3.24) 
/ dn ■ 


'on 


x)] 4/(x) — x(x)^D4/(x) 


= 0 , 


and one obtains 


(JDx{x)^ X (j4/(x) — x(x)^(j X D4/(x) = V x (x(x)^a4/(x)) . (3.31) 

Applying Stoke’s theorem and using the fact that the boundary of a boundary is an 
empty set (i.e. d{dfl) = 0) we arrive at 


dn ■ 


'an 


/ - , .\t 


yDx{x)j X (j\k(x) — x(x)V X D'h(x) 


Ian 


dn-V X (x(x)'^a'h(x)) = / dl ■ x(x)'''a'h(x) = 0. (3.32) 


ia{dn) 


The Hermiticity condition eq.(|3.30|) can thus be expressed as 
/ dn ■ 


'an 


(j5x{x) — ia X Dx(x)^ ^(x) — x(x)^ ^D4/(x) — ia x DT(x) j 


= 0 . 
(3.33) 
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Using the self-adjointness condition eq.(3.27), we then obtain 


/ cPx ~ = 0. (3.34) 

JdQ 

Here we have used the fact that 7 ( 0 ;) is Hermitean. 

The matrix 'j{x) results from the matrix A(T) in the non-relativistic limit. Rewrit¬ 
ing the self-adjointness condition eq.(3.27) as 

7 (T)T(T) -|- n{x) ■ aa ■ D'^{x) = 'y{x)'^{x) + 2Mcin{x) ■ aX{x)'^{x) = 0, (3.35) 


one arrives at 

7 (T) = —2Mcin{x) ■ aX{x). (3.36) 

Using the anti-Hermiticity of n{x) ■ aX(x), one concludes that 'y(x) is indeed Her¬ 
mitean. 


4 The Hydrogen Atom in a Spherical Cavity 


In this section we place a hydrogen atom at the center of a spherical cavity and 
investigate the effect of the boundary conditions on the accidental symmetry. Again, 
we hrst study the Schrodinger equation before we proceed to the Dirac and Pauli 
equations. 


4.1 The Non-relativistic Schrodinger Atom in a Spherical 
Cavity 


Let us consider the Schrodinger Hamiltonian for the hydrogen atom, 


1 

H = -A-= 

2M r 


2M 


dl + -a, - — - 


For the wave function we make the factorization ansatz 

T(f) = 

which leads to the radial equation 


(dr + Idr - ^ V'nz(r) = i7V'nz(r). 


When we write the energy as 


E = 


Me^ 


(4.1) 


(4.2) 


(4.3) 


(4.4) 
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for the bound state spectrum in the inhnite volume n is quantized in integer units. 
Inside a spherical cavity, on the other hand, n can take arbitrary values, and one 
obtains 

{^) {^) “p 

where /na) is an associated Laguerre function. Imposing the most general 

boundary condition, 7 (a;)\I'(a;) + n{x) ■ V\I'(a;) =0, one obtains 

'yi>ni{R) + dri>ni{R) = 0 , (4.6) 


where R is the radius of the spherical cavity. Consequently, the quantized energies 
result from the equation 


'jna 

~Y~ 


I l^\ 21+1 
2^ 2r) 


_ t21+2 
^n-l-2 



0 . 


(4.7) 


The resulting hnite volume spectra for 7 = oo (i.e. Dirichlet boundary conditions) 
and for 7 = 0 (i.e. Neumann boundary conditions) are illustrated in Figs. [^and|^ 
respectively. 



Figure 1: Spectrum of the Schrodinger hydrogen atom centered in a spherical cavity 
with the standard Dirichlet boundary condition (i.e. 7 = ooj as a function of a/R. 
The energies of states with guantum numbers n = 1,2,3 are given in units of Me‘^. 
The dotted lines represent the spectrum of the infinite system (cf. f23f). 


As we discussed in Section 2 . 1 , in the inhnite volume the non-relativistic hy¬ 
drogen atom enjoys an accidental S'0(4) symmetry gener ated by the angular mo¬ 
mentum L together with the Runge-Lenz vector R of eq.(2.2), which results in an 


14 

















Figure 2: Spectrum of the Schrodinger hydrogen atom centered in a spherical cavity 
with Neumann boundary condition (i.e. 'y = 0) as a function ofa/R. The energies 
of states with guantum numbers n = 1,2,3 are given in units of Me^^. The dotted 
lines represent the spectrum of the infinite system (c.f. 


n^-fold degeneracy of the bound state spectrum. The spectrum of hydrogen conhned 
to a spherical cavity, on the other hand, no longer shows the accidental degeneracy. 
Since the boundary condition does not violate rotation invariance, L is still con¬ 
served, but R is not. This is because the application of the Runge-Lenz vector leads 
outside the domain of the Hamiltonian |23] , as was already pointed out in [25l [27] 
for the Dirichlet boundary condition with 7 = 00 . Relying on rotation invariance, 
we restrict ourselves to states T(r) = , (p) with the maximal m = 1. On 

such states the raising operator = R^ + iRy acts as [23] 


R+T(r) = 


/ + ! 
M 


dr'tpniir) + e 


/(/ + !) 

Mr 


f’niir) 


Y, 




{d,T) 


= Xn,i+i{r)Yi+i^i+i{e,p). 


(4.8) 


After an application of R the wave function remains in the domain of the Hamilto¬ 
nian only if the new wave function y;„ ;_|_i(r) also obeys the boundary condition 


dXrijZ-l-l (R) T drXn,l+l{R) 0. 


(4.9) 


Inserting eq.(4.8) into this relation and using the boundary condition x'f’niiR) + 
dr'ipniiR) = 0 as well as the radial Schrodinger equation (4.3) for fjniir), one hnds 


XXn,l+l{,Rf T 9rXn,l+l{,Rf 

[±1 + 
M [ ^ V R/ R^ 


2Me^ 

R 


-2ME 


'IpnliR)- 


(4.10) 
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Since the right-hand side of eq.(4.10) does not vanish independent of the energy E, 
the new wave function Xn,/+i(’") that results from the application of R+ on 
does not obey the boundary condition eq.(4.9), and hence lies outside the domain of 
the Hamiltonian. Consequently, the Runge-Lenz vector R leads out of this domain, 
and thus no longer represents an accidental symmetry. 



a/R 



arctan(Y R) 


Figure 3: Top: Energy of s- and d-states for the Schrodinger hydrogen atom in a 
spherical cavity with 7 = 1/a as a function of a/R. There is an accidental degeneracy 
for R = 2a. Bottom: Energy of s- and d-states for a hydrogen atom in a spherical 
cavity with radius R = 2a as a function o/arctan( 7 R). The energies of the 3s and fd 
states are very similar for all values of but are identical only for 7 = 2 /R = 1 /a 
and for 7 = ±cx). The energy is given in units of Me^. 
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Interestingly, for 7 = cx) a remnant of the accidental symmetry persists even 
in a finite volume, however, only for the special cavity radius R = (Z + 1)(Z + 2)a 
PSI ITT] , Then a state with angular momentum I is degenerate with a state of 
angular momentum 1 + 2. In particular, for R = 2a, the states 2s and 3d, 3s and 
4d, 4s and 5d, etc. form multiplets of 1 + 5 = 6 degenerate states, while for R = 6 a, 
the states 3p and 4f, 4p and 5f, 5p and 6 f, etc. form multiplets of 3 + 7 = 10 
degenerate states. These accidental degeneracies are again due to the Runge-Lenz 
vector. They arise because, for the special values R = (/ + !)(/ + 2)a, the operator 
R^ maps 'ipnii'f) back into the domain of the Hamiltonian. Indeed, one can show 
that R\i)ni{r)Yii{e,(p) = Xn,i+ 2 {r)Yi+ 2 ,i+ 2 { 0 ,+) with 


Xn,i+ 2 {r) 


2/ + 3 / 2 (Z + 1)(Z + 2) 
M \ Mr 

■(Z + l)(Z + 2 ) /Z(2Z + 3) 

M V 


dri^niir) 




Z(Z + 1) \ 
Mr J 


(4.11) 


Since for 7 = cx) the wave function obeys the Dirichlet boundary condition 'ipni{R) = 
0 , for 


2 _ (Z + 1 )(Z + 2 ) 
MR 


0 i? 


(Z + l)(Z + 2) 
Me2 


(Z + 1) (Z + 2)a, 


(4.12) 


the wave function Xn,i+ 2 {f') indeed obeys the same condition. Remarkably, the same 
accidental degeneracy arises for R = (Z + 1)(Z + 2)a and 7 = 2/i? |26]. This fol¬ 
lows from the fact that then x'^nii.R) + dr'^nii.R) = 0 indeed implies 77 ^ ,1+2{R) + 
drXn, i+ 2 {R) = 0. The remnant accidental degeneracies of the non-relativistic hydro¬ 
gen atom in a spherical cavity are illustrated in figure 


4.2 The Relativistic Dirac Atom in a Spherical Cavity 


Let us now consider the Dirac equation in a spherical cavity, i.e. we impose the 
boundary condition 


( ^3(^) 

V ^4(^) 



x\ = R, 


(4.13) 


where a ■ erX{x) is an anti-Hermitean matrix. Again, we want to maintain the 
spherical symmetry, [H, J] = 0. We thus act with on a general Dirac spinor \h(a;) 
and demand that = J“\I'(x) still obeys the boundary condition 


A(f) 


T2(f)' 




Ti(f) 

T2(f) 


\ 

^4ixy J 


( ^3(^) 

V ^ 4 (^) 


r\{x) 


Ti(T) 

^2{x) 



(4.14) 
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This hence implies that [J“, A(a;)] = 0. Similarly, the boundary condition should 
also maintain the parity symmetry P, such that \h(a;)' = P'^{x) should again obey 
the boundary condition 


A(f) 


^ 2 (x)' J 


\{x) 


Ti(-f) 

T 2 (-T) 


V -d'4(-T) J 


V ^4(x)' ) 


-A(-x) 


Ti(-f) A 


(4.15) 


This implies the restriction A(a;) = — A(—T). Rotation invariance and parity limit 
the anti-Hermitean matrix that characterizes the boundary condition to 


A(a;) = iua ■ Cr, 


(4.16) 


such that a ■ er^{x) = izzl is indeed anti-Hermitean. Let us further investigate 
whether K remains a symmetry in the hnite volume. For this purpose we act 
with K on the Dirac spinor 4/(T), and we ask once more whether the new spinor 
^/(T)' = K'^{x) still obeys the boundary condition. This is the case only if 


A(f) 


d^i(x)' A 

'^2ixy J 




which is indeed satished since 


^3(x)' A 
d^4(f)' J 


(4.17) 


{A(a;), a ■ L + 1} = {iua ■ e)., a ■ L -|- 1} = 0. (4-18) 


Next, we want to investigate whether the Johnson-Lippmann operator also leaves 
the boundary condition invariant. Hence, we now ask whether 4/(1;)' = H4/(l;) also 
obeys the boundary condition, which is the case only if 


A(f) 


4/i(l;)' 

4/2(T)' 


Using [A(T), a ■ e).] 
identity holds 


X{x)^^{a ■ L + l)a ■ p — X{x)aa ■ e7 


A(f) 


la 

McR 


T 3 (T)' 

4/4(T)' 


{a ■ L + l)A(T) 


'Fi(T) 

'F 2 (T) 


la 

McR 


(jj ■ L P P) 


-^^(d ■ L -|- l)a ■ pX{x) — aa ■ erX{x) 


'Fi(T) A 
'F 2 (T) ) ' 

(4.19) 


0 , this condition would indeed be satished if the following 


{A(T), (d - L + l)d-p} = ^ 
iv{a ■ er,{a ■ L + l)d ■ p} = 


(^X{x){a ■ L + 1)A(1;) + a ■ L + 1^ 

|(z 22 + i)(d-r+i). 


(4.20) 
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It is straightforward to show that 


i{a ■ Cr, {a ■ L + l)a ■ p} = {a ■ L + \)i[a ■ p, a • e^] = —{a-L + 1)^, (4.21) 

R 

which implies that A maintains the boundary condition if 

= a-L + l = k = ±(^j + ^. (4.22) 

Obviously, this relation cannot hold as an operator identity. Still, it can be satis- 
hed for states with an appropriate /c-value and for specihc values of u. One might 
then conclude that, like in the non-relativistic case, a remnant accidental degener¬ 
acy arises for particular states and for specihc values of the self-adjoint extension 
parameter. However, the situation is more subtle. In particular, since {K, A] = 0, 
an application of A changes the sign of k. This means that, for a hxed value of 
z/, eq.(4.22) is no longer satished for —k, and thus another application of A will 


inevitably lead out of the domain of H. Since repeated applications of H or K, on 


the other hand, leave the wave function inside the domain, eq.(2.16), which relates 


to H'^ and K"^ in the inhnite volume, is no longer satished in a hnite volume. In 
fact, H and A no longer commute as operators restricted to the domain of H. This 
implies that A no longer generates an accidental symmetry, and the corresponding 
degeneracy is lifted in a hnite volume. 

Let us now investigate the spectrum of relativistic hydrogen in a spherical cavity 
in more detail. First we make a separation ansatz for the wave function 

/ 4'i(f) \ 

T2(f) ^ 

4 ' 3 (^) 

V 4'4(x) 

For k = j + ^ we have Ia = j — \ and Ib = i + while for k = —(j -f 4), = j + 4 

and Ib = j — For j = I + \ the spin-angular functions are given by 


T(T) = 


(4.23) 




^ + is + L, fn f 1 ^ “43 + i 


2,+i 0 ' + 


2 / + 1 




0 


while for j = I — \ 


(4.24) 






2,+ l 1 


(4.25) 

Here Yim{9,ip) are the usual spherical harmonics. The radial equations then take 
the form 

.2X / 


r 


Me - i/jA{r) - c[dr + 


r 


^/’B(r) = EtjjAir), 


-Mc^ - i’Bir) + c (dr + ) i’Air) = EipBir). (4.26) 
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The cavity boundary condition is given by 

iucr ■ a'ii)A{R)yjj3iA^, = -w'ipA{R)yjhiB{o, = i'ipB{R)yjj3iB{o, => 
v^Pa{R) = -^b{R). (4.27) 

The resulting hnite volume spectra for v = oo (i.e. Dirichlet boundary conditions) 
and for i/ = 0 (i.e. Neumann boundary conditions) are illustrated in Figs. |^and|^ 
respectively. In order to make the effects easily visible we have chosen very large 
unphysical values of a. In both cases, we see that the states 2Si/2 and 2Pi/2, which 
are accidentally degenerate in the inhnite volume, are split in the hnite cavity. With 
Dirichlet boundary conditions (z/ = cx)) the energies increase with decreasing cavity 
radius, while for Neumann boundary conditions [v = 0) the energy of some states 
decreases. 


4.3 The Non-relativistic Pauli Atom in a Spherical Cavity 

Finally, we proceed to the non-relativistic limit of the Dirac equation, which leads 
us to the Pauli equation 


H = My + 


pc + eA{x) 
2 Mc 2 


— e<h(T) + pa ■ B{x). 


(4.28) 


Since we have no external magnetic held, at least to the leading order we are working 
at, spin decouples and the Pauli equation reduces to two copies of the Schrodinger 
equation. In particular, we have neglected the sub-leading spin-orbit couplings in 
the Hamiltonian. Still, as we discussed in Subsection 3.3, the spin explicitly enters 
the boundary condition 


7(x)\['(T) -f n(T) ■ V\k(T) — Z(j X VT(x) 


= 0 , 


(4.29) 


which thus leads to a spin-orbit coupling induced by the cavity wall. The Her- 
mitean matrix 7 (T), which results from the anti-Hermitean matrix A(a;) in the non- 
relativistic limit, is given by 


7 (T) = —2Mcin{x) ■ aX{x) = 2Mch'l. 


(4.30) 


For j = I + h the boundary condition of eq.(4.29) reduces to 


while for j = I — \ 


'fijnliR) + dri^nliR) “ —ll^nl{R) = 0 , 

it 


'yi>nl{R) + dri>nl{R) + ^-^^i>nl{R) = 0 , 


(4.31) 


(4.32) 
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Figure 4: Spectrum of the Dirac hydrogen atom centered in a spherical cavity with 
the standard Dirichlet boundary condition (i.e. u = oo) for a = 1 (top) and for 
a = (bottom), as a function ofa/R. The energies of states with guantum 

numbers n = 1 and 2 are given in units of M. The dotted lines represent the 
spectrum of the infinite system. 

The spectrum of the Pauli equation for Neumann boundary conditions (i.e. 7 = 0) is 
illustrated in Fig. It differs from the spectrum of Fig. [^because now the angular 
momentum enters the boundary condition. 

Let us now address the question of remnant accidental degeneracies that persist 
even in a hnite volume. As we discussed in Subsection 4.1, for the Schrodinger 
atom these exist for R = (/ + !)(/ + 2)a and 7 = 2/R or 00 , because, in these 
cases, two applications of the Runge-Lenz vector lead back into the domain of the 
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Figure 5: Spectrum of the Dirac hydrogen atom centered in a spherical cavity with 
Neumann boundary condition (i.e. v = Oj for a = 1 (top) and for a = 

(bottom), as a function ofa/R. The energies of states with guantum numbers n = l 
and 2 are given in units of M. The dotted lines represent the spectrum of the infinite 
system. 


Hamiltonian. For the Dirac atom in a finite volume, on the other hand, repeated 
applications of the Johnson-Lippmann operator do not lead back into the domain 
of the Hamiltonian, and thus the accidental degeneracy is lifted in a finite volume 
(c.f. Subsection 4.2). 


As we will see now, just as for the Schrodinger atom, for the Pauli atom some 
accidental degeneracies persist in a finite cavity. We start again from eq.(4.11), take 
its derivative with respect to r, and then eliminate by using the radial 
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Figure 6 : Spectrum of the Pauli hydrogen atom centered in a spherical cavity with 
Neumann boundary condition (i.e. 7 = Oj as a function of a/R. The energies of 
states with guantum numbers n = 1,2,3 are given in units of Me*. The dotted lines 
represent the spectrum of the infinite system. 


Schrodinger equation (4.3). We then impose the boundary condition eq.(4.31) when 
j = I + or eq.(4.32[) when j = / — 4. Finally, we ask under what circumstances the 


wave function Xn,i+ 2 {r) of eq.(4.11) satishes the corresponding boundary condition 


XXn,l+2iR) + drXn,l+2iR) 

when j' = / + 2 + 4 = / + |, or 


1 + 2 


Xn,l+2{R) = 0, 


(4.33) 


7 Xn,«+ 2 (i?) + drXn,l+2{.R) + ^^^Xn,Z +2 (i?) = 0, (4.34) 

when j' = / + 2 — 4 = / + |, for all values of the energy E. For j = I + \ and 
i' = I + \ this turns out to be the case for 


(/ + l)(/ + 2)(2/ + 5)( 
2/+ 3 


7 = 


■ 2 (/+ !)(/ + 2 )a’ (/ + l)a' 




For i = I + \ and j' = / + |, on the other hand, the conditions can be satished only 
when = 0 and 7 = 0 . For j = I — \ and f = I + \ one obtains 


R = 2{l + l){l + 2)a, 


1 

2 (/ + l)(/ + 2 )a 



3 _ 3 

2 (/ + l)(/ + 2 )a “ R' 
(4.36) 
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1 


(4.37) 


and, finally, for j = / — ^ and j' = I + \ one finds 

{I -\- 1)(^ + 2)(2/ + l)a 1 

“ 2/+ 3 ’ ~2(/+ !)(/ +2)a’ ^ ~ + 2)a' 

The remnant accidental symmetries for the Panli hydrogen atom are illnstrated in 
Fig. [71 









R = 12a _ P-states 

y=-l/12» .,=1/40 0 
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r ifui 





arctan(Y R) 
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-- 
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Fignre 7: Top left: Energy of P- and F-states with n = A and 5 for the Pauli 
hydrogen atom in a spherical cavity with 7 = — ]^ as a function of a/R. There are 
accidental degeneracies (indicated by the vertical dashed lines) of the 4 P 3/2 with the 
4^7/2 for R = y a, as well as of the 4 P 1/2 with the 4 P 7/2 and of the 5Pi/2 with the 
5 F 7/2 states for R = 12a. Top right and bottom: Energy of P- and F-states for the 
Pauli hydrogen atom in a spherical cavity as a function o/arctan( 7 i?) at fixed cavity 
radius: R = (top right), R = 12a (bottom left), and R = (bottom right): 
The energies of the high-lying P- and F-states are very similar for all values of 7 , 
but are identical only for the special values indicated by the vertical dashed lines. 
The energy is given in units of Me'^. 


^We like to thank D. Banerjee for valuable help with deriving eqs.(4.35), (4.36), and (4.37). 
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5 Conclusions 


We have investigated the accidental symmetry in the relativistic hydrogen atom 
conhned to a spherical cavity with perfectly reflecting boundary conditions. In the 
inhnite volume the Johnson-Lippman operator A, which is the relativistic analog 
of the Runge-Lenz vector R, commutes with the Hamiltonian, thus leading to acci¬ 
dental degeneracies in the energy spectrum. When the system is placed in a finite 
volume, the accidental symmetry is lifted. This is because a repeated application of 
the Johnson-Lippman operator leads out of the domain of the Hamiltonian. In the 
non-relativistic case, for specific values of the cavity radius R and the self-adjoint 
extension parameter 7 , a repeated application of R leads back into the domain of 
the Hamiltonian, and thus some accidental degeneracy persists even in a hnite vol¬ 
ume. Because repeated applications of A do not lead back into the domain of the 
Hamiltonian, this is not the case for the Dirac equation. Interestingly, in the Pauli 
equation with spin, which results as the non-relativistic limit of the Dirac equation, 
the boundary conditions induce non-trivial spin-orbit couplings. Remarkably, in 
this case, for specific values of the cavity radius R and the self-adjoint extension 
parameter 7 , a repeated application of R again leads back into the domain of the 
Hamiltonian, and thus some accidental degeneracy persists in a finite volume. 

Our investigation shows that the subtle issues of Hermiticity versus self-adjoint- 
ness and the domain structure of quantum mechanical Hamiltonians hold the key 
to understanding in detail what happens to a hydrogen atom when it is confined in 
a spherical cavity. The same is true for other accidental symmetries, for example, 
for the Landau level problem [S], for a particle on a cone |23], or for the harmonic 
oscillator [221121] • This underscores that the theory of self-adjoint extensions is not 
just a mathematical curiosity, but of great importance for understanding the physics 
of conhned systems. 
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